Avoided Band Crossings: 
Tuning Metal-Insulator Transitions in Chaotic Systems 
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We show that avoided crossings of energy bands may give rise to a variety of phenomena such as 
transitions from metal to insulator and vice versa, changes in localization lengths, and changes in 
the fractal dimension of energy spectra. We explain the occurrence of these phenomena in the kicked 
Harper model under classically chaotic conditions and predict them to occur in other systems. 
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Spatially periodic quantum systems give rise to energy 
bands and extended eigenfunctions, according to Bloch's 
theorem. If this periodicity is broken by a weak pertur- 
bation, e.g., by disorder or a (magnetic) superlattice, the 
spectrum remains band-like. Typically, one neglects the 
coupling of these "bands" and studies them in a one- 
band approximation. Prominent examples from solid 
state physics for such one-band models are the Anderson 
model for disordered systems pi] and the Harper model 
for Bloch electrons in a magnetic field |^,|| , which display 
such interesting properties as localized eigenfunctions or 
metal- insulator transitions. As a function of an exter- 
nal parameter some of these "bands" will move up in 
energy whereas others move down, thus giving rise to 
band crossings. At these crossings the coupling of bands 
can no longer be neglected and will lead to avoided band 
crossings (ABCs) in analogy to the well-studied avoided 
level crossings. The question then arises: Will these cou- 
plings of bands and in particular the ABCs change the 
properties derived in the simplifying one-band models? 

A seemingly unrelated problem has occurred in the 
quantum chaos literature, while studying the influence 
of classical chaos on quasiperiodic quantum systems. 
Specifically, in the kicked Harper model (KHM) [[| pl| . 
it was found numerically that in the regime, where the 
classical limit is chaotic, extended eigenfunctions show up 
for parameters where exclusively localized ones were ex- 
pected and vice versa [|]||. In addition, chaos was found 
to change the fractal dimension of the spectrum ||. In 
the absence of an explanation these puzzling phenomena 
have remained mysterious, and one may wonder: How 
can classical chaos induce such drastic changes in these 
quantum properties? 

In fact, it will turn out here that the main impact of 
classical chaos is to generate ABCs and thus this ques- 
tion for the quantum chaos problem will be reduced to the 
previous question posed for the solid state problems. In 
the main part of this paper we therefore analyze avoided 
crossings of one-band models. Although simple gener- 
alizations of avoided level crossings, these ABCs unex- 
pectedly generate transitions from localized to extended 



eigenfunctions and vice versa and changes in the local- 
ization lengths of eigenfunctions as well as in the frac- 
tal dimensions of spectra. We show that ABCs cause 
these phenomena by effectively changing the parameters 
of the one-band models involved. We will demonstrate 
that these general findings can be applied to the KHM 
and thereby give the first explanation for the phenomena 
that were observed numerically in this model. Applica- 
tion to the old problem of Bloch electrons in a magnetic 
field Q leads to the prediction of dramatic effects, e.g., a 
wave packet spreading in one direction will, due to ABCs, 
spread in the perpendicular direction. These effects give 
rise to experimental consequences for transport measure- 
ments in semiconductor nanostructures and will be de- 



scribed in detail elsewhere (16). In the future we expect 
applications to coupled one-dimensional chains and many 
other systems. 

Let us start by studying an ABC when there is no 
perturbation and the system is still periodic. Then each 




FIG. 1. An ABC with twist modeled by the Hamilto- 
nian ([j]), with E' = — cosfc, E" — cosfc, and e = 0.6. The 
eigenvalues form two energy bands for the set of all Bloch 
phases k. Arrows indicate in which direction the energies 
change as k is varied from to n. Solid and dashed lines rep- 
resent the eigenvalues for k — and k = ir, resp. The inset 
shows an ABC without twist for the case E' = E" = — cos k 
with e = 1.4. 
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band is described by a dispersion E(k) depending on the 
Bloch phase k 6 [— ir, w]. Now, an isolated avoided cross- 
ing of two bands can be modeled by the Hamiltonian 

H -[ e E»(k)-\)' (1) 

where e is the coupling parameter. For any fixed Bloch 
phase fc, H simply models an avoided level crossing. 
When looking at the entire band there are two cases of 
interest: If the dispersion relations of the two bands have 
the same general dependence on k, namely 4§ > for 
k G [0,7r], the resulting spectrum resembles that of a 
broadened avoided level crossing (see inset of Fig. 1). In 
the alternative case of opposite general dependencies on 
k, namely > and < for k <E [0, n], each band 
is twisted (Fig. 1). 

Two subsequent ABCs with twist (see Fig. 2a) can be 
modeled by the Hamiltonian 

/ E+(k)+A-X e \ 

H=\ e E°{k) + X e 

\ e E-{k)-A-X j 

(2) 

where a rising band E subsequently crosses two decreas- 
ing bands E~ and E + with offsets ±A. For small cou- 
pling e the two twists of the intermediate band are well 
separated, whereas for increasing e they first come closer 
whereby the band flattens, which will be important later 
on, and finally the twists annihilate (Fig. 2). 

We now show the dramatic effect of ABCs in the pres- 
ence of a weak perturbation. This we model by the 
Hamiltonian of Eq. (||) with the bands E ±,0 (k) replaced 
by Hamiltonians iJ ±! °, leading to the Hamiltonian 

/H+ + A-X e \ 

H = [ e H° + X e . (3) 

\ e H- -A-X J 

For simplicity, we will focus on tight-binding Hamiltoni- 
ans H±>° = £ B V*=-°at °n+E„ t±>° (4 +1 a„ + a]^) 




FIG. 2. Two subsequent ABCs with twist generated by the 
Hamiltonian (|^) with E°(k) — — cosfc, E ± (k) = cosfc, and 
offset A = 3 for increasing coupling e. For small e the twists 
are well separated in A (a). Upon increasing e the two twists 
of the intermediate band approach each other and reduce its 
bandwidth (b,c). If e is increased further, the twists annihilate 
(d). 
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FIG. 3. The spectrum of Eq. (|3j) where H represents 
the Harper model in the metallic regime (V°/t° = 1.2, 
a = 5/8) and H ± has zero on-site potential and hopping 
terms t — —t° leading to twists. The coupling parameter 
e = 1.5 corresponds to Fig. 2c. At the bottom typical eigen- 
functions are shown where one clearly sees a transition from 
extended eigenfunctions in the uncoupled regime (|A| 2> 1) 
to localized eigenfunctions in the regime between the twists 
(A = 0). 

on one-dimensional lattices with spatially varying on-site 
potential V^' and with hopping terms i ±! °. For zero 
on-site potentials (V±'° = 0) Eq. © reduces to Eq. (§) 
studied in Fig. 2. 

As a first numerical example we assume H° accord- 
ing to Harper's model for Bloch electrons in a magnetic 
field H in the metallic regime, i.e., we take the on-site po- 
tential V° = V° cos(27rcm) and V°/t° < 2, where a is the 
magnetic flux per unit cell. For H we assume vanishing 
on-site potential = and, in order to generate twists, 
hopping terms t = —t°. Fig. 3 corresponds to Fig. 2c, 
where the three tight-binding bands now show fine struc- 
ture, while on a coarse scale they display the same twists 
and flattening as in Fig. 2c. For large |A| the individual 
spectra of H ± and H° are uncoupled and all eigenfunc- 
tions are extended. Surprisingly, in the A region between 
the twists the intermediate band has localized eigenfunc- 
tions only. If for H° we assume the Anderson model for 
disordered systems 0, where V® is a uniform random 
variable in the interval [-V°/2, V°/2] , for A = we find 
a decreasing localization length of eigenfunctions as the 
coupling parameter e is increased (Fig. 4a, squares). As 
a third example for H° we assume the Fibonacci chain 
model |r7f| , a one-dimensional model for quasicrystals 
with potential terms taking on the values V° and — V° ac- 
cording to a Fibonacci sequence. Remarkably, at A = 
increasing e lowers the fractal dimension of the energy 
spectrum (Fig. 4b, squares). These three systems show 
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FIG. 4. a) Average inverse participation ratio 
£ = Yin IV^n I* °f the eigenfunctions tp of the intermediate 
band of Eq. ([]) for X = 0, A = 10, and different values of e. 
Here 7/° represents the Anderson model (V° = 5, t° = 1) and 
has zero on-site potential and hopping terms t = — 1 
(squares) and t — 1 (diamonds), resp. Solid (t — —1) and 
dashed (t = 1) lines are the results for an isolated Ander- 
son model with the effective terms V® and i of perturbation 
theory (Eqs. (^) and (^|)). An increasing (decreasing) inverse 
participation ratio corresponds to a decreasing (increasing) 
localization length. 

b) Same as (a) showing the fractal (box-counting) dimension 
for H° representing the Fibonacci model (V° = 1.5, t° = 1, 
A = 100). 

the opposite of what one would expect naively: The 
coupling of a given tight-binding model (H°) to mod- 
els having extended eigenfunctions (H ) in these three 
cases leads to localization, reduced localization lengths, 
and reduced fractal dimensions, respectively. How is that 
possible? 

One can try to understand these effects of ABCs in- 
tuitively: The reduced width of the intermediate band 
in Fig. 2c resulting from subsequent ABCs with twists 
corresponds to reduced hopping terms in a tight-binding 
Hamiltonian [ fL8| . More generally, one may expect that 
the ratio of potential terms to hopping terms is increased 
by these ABCs with twist. This would naturally explain 
the observed results, as for increasing V°/t° the Harper 
model makes a transition from metal to insulator Q , the 
Anderson model has a decreasing localization length jp , 
and the energy spectrum of the Fibonacci model has a 
decreasing fractal dimension |l9(| . 

Further insight can be obtained from a perturbative 
treatment for A = 0. By a similarity transformation we 
obtain H = S~ 1 HS, where 



S = 



[l-e 2 /2A 2 ) 
e/A 
e 2 /2A 2 



-e/A 
(1 - e 2 /A 2 ) 
e/A 



e 2 /2A 2 
-e/A 
<l-e 2 /2A 2 ) 



(4) 



This choice of S leads to renormalized Hamiltonians H ±: ° 
analogous to Eq. (|^) that are decoupled up to order e 2 / A 2 
for \V±'°\, |i ±,0 | < A. The Hamiltonian H° describing 
the intermediate band is still tridiagonal with on-site po- 
tential V® and hopping terms i° up to order e 2 / A 2 given 
by 



i° = t°(l-4e 2 /A 2 ) 



(5) 



For the systems studied above (where = V~ = 0) we 
find that the ratio V®/t° is effectively increased by the 
factor (l + 2e 2 /A 2 ). This substantiates the intuitive ex- 
planation of the preceding paragraph. E.g., in the Harper 
model ABCs with twist may induce a transition from ex- 
tended states (V°/t° < 2) to localized states (V°/t° > 2). 
For two subsequent ABCs without twist, i.e. when the 
hopping terms ^ and t° have the same sign, the above 
perturbative treatment yields 



V° = V°(l-2e 2 /A 2 ) + (V+ + V-)e 2 /A 2 



t v = t" 



(6) 



instead. Thus, without twists the ratio V®/t° effectively 
decreases by the factor (1 — 2e 2 /A 2 ). In this case we 
therefore expect the opposite effects on eigenfunctions 
and spectrum. Using the renormalized terms of Eqs. (|5|) 
and (|6|) we can compute the localization length for an 
isolated Anderson model and the fractal dimension Do 
for an isolated Fibonacci band. Comparison with the nu- 
merical results in Fig. 4a) and b) gives a good agreement 
and thus confirms our explanation in terms of renormal- 
ized tight-binding parameters. In contrast to the case 
of ABCs with twist, here an intuitive understanding is 
missing, and one has to rely on the above perturbative 
treatment. The phenomena associated with ABCs ana- 
lyzed above also do occur when there is just a single ABC, 
although in a more complex way. Thus they should be 
generic for all models with ABCs. 

Now we will apply these general results to the well- 
studied kicked Harper model (KHM) @-|lE| . This model 
combines the quantum mechanical properties of Harper's 
equation with an onset of chaos in classical phase space. 
It shows similar phenomena , as those discussed for 

ABCs, but they have remained unexplained so far. The 
KHM is given by the Hamiltonian 



H = L cosp - 



K cos x 8 (t 



(7) 



with p — —ifid/dx and an effective Ti. For small K and L 
its spectrum is close to the Harper spectrum for a = %/2ir 
and V°/t° — 2K/L up to a scaling by a factor propor- 
tional to K H (Fig. 5). In particular, the spectrum clus- 
ters into subbands each of which may be described by 
a Harper-like model [EQ]. Upon increasing the parame- 
ters K and L these subbands make band crossings. As 
soon as the classical limit is non-integrable (i.e., mixed or 
fully chaotic) such crossings are avoided. Therefore our 
general findings on ABCs apply and one may expect all 
the phenomena discussed above to occur. Indeed, this 
is the case. For example, for the parameters of Fig. 5 
transitions from extended to localized eigenfunctions can 
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FIG. 5. Transitions of the quasienergy spectrum of the 
kicked Harper model for increasing kicking strength K and 
L/K — 4/7 as the classical phase space becomes increasingly 
chaotic (insets on top). Here h/2ir = 8/61 is a rational ap- 
proximant of 1/(6 + a a) where oq = 1.618... is the golden 
mean and the Bloch phase in p-direction is varied. While for 
small K and L all eigenfunctions are extended as expected (a 
typical one is shown at the bottom), for a chaotic phase space 
and stronger kicking K one simultaneously finds localized and 
extended states (bottom) depending on whether or not ABCs 
and twists have led to flattened bands. 

be observed as were numerically found in Refs. HI). 
They occur along with ABCs, however, in a more com- 
plex way than in the simple model system of Eq. (|^). 
Tuning these metal-insulator transitions necessarily leads 
to states right at the transition point, which, e.g., in 
the Harper model correspond to a singular continuous 
spectrum. We therefore also understand the findings of 
Borgonovi and Shepelyansky Jll|, that for the dual pa- 
rameters of Fig. 5 at K = 4 and L = 7 part of the 
spectrum seems to be singular continuous Simi- 
larly, the increase of the fractal dimension of the spec- 
trum for increasing K — L @] can be understood on the 
basis of ABCs as analyzed above in the Fibonacci model 
(Fig. 4b). 
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